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ABSTRACT 

Weak gravitational lensing provides a potentially powerful method for the detection of clus- 
ters. In addition to cluster candidates, a large number of objects with possibly no optical 
or X-ray component have been detected in shear-selected samples. Determining the nature 
of these so-called "dark" lenses is an important step towards understanding the reliabil- 
ity of shear-selection te chniques. We develop an analytic model to investigate the claim of 
I Weinberg & Kamionkowski (2002) that unvirialised protoclusters account for a significant 
number of dark lenses. In our model, a protocluster consists of a small virialised region sur- 
rounded by in-falling matter. We use a simple model for the density profile that assumes the 
Navarro-Frenk- White form inside of the virial radius and a power law p ~ r~" outside. We 
find that, in order for a protocluster to simultaneously escape X-ray detection and create a 
detectable weak lensing signal, it must have a small virial mass 10^'^ M©) and large total 
mass (~ 10^^ Mq), with a relatively flat density profile outside of the virial radius (a ^ — 1). 
Such objects would be characterized by rising tangential shear profiles well beyond the virial 
radius. We use a semi-analytic approach based on the excursion set formalism to estimate the 
abundance of lensing protoclusters with a low probability of X-ray detection. We find that they 
are extremely rare, accounting for less than 0.4 per cent of the total lenses in a survey with 
background galaxy density n — arcmin^^ and intrinsic ellipticity dispersion — 0.3. 
Their abundance decreases significantly if flat density profiles outside of the virial radius are 
not common. We conclude that lensing protoclusters with undetectable X-Ray luminosities 
are too rare to account for a significant number of dark lenses. 

Key words: cosmology: theory - large-scale structure of the universe 
; 1 INTRODUCTION 

' The abundance of collapsed dark matter haloes in the Universe yields a wealth of information on the dynamics of structure formation. The 
number density of galaxy clusters a nd its time-evolution can be u s ed to probe the normalization of the linear power spectrum, ag,, and the 
density parameters and (e.g..'Lilie"l992VWhite et al.'l993';'Cen & Ostriker'l 9"94| ; lEke et al.ll996l : lHenrvll997l : lBahcall & Fanll998l : 
[viana & Liddle 1999; Holder et al. 2001; Reiprich & Bohringer 2002; Dahle 2006). In addition, since the linear growth of overdensities and 
comoving volume-element are sen sitive to the dark energy e quation-of-state parameter w = P/ p, the cluster abundance can be used to 
constrain dark energy (e.g.. iHaiman et alJl200ll ; lMaiumdar & Mohj|2004lMantz et alj|2008l) . In order to effectively use the cluster number 
counts as a cosmological tool, one must be able to create a complete catalog out to high red shifts. To da te, most studies aiming to complete 
this task have relied on X-ray based selection and constraints on cluster masses (see however [Dahlel20 0^. Mass estimates derived from their 
X-ray emission require the additional assumption of hydrostatic equilibrium for the gas, which is not robust. 

Galaxy clusters are the most recently assembled structures in the Universe. One should therefore expect to find a large number of them 
in a dynamically unrelaxed state. In recent years, much effort has been devoted to developing techniques that utilize weak-lensing in cluster 
surveys. Gravitational lensing offers a method to measure masses that is independent of their dynamical state. In addition to improving the 
accuracy of mass measurements, the weak distortion in the shapes of background galaxies may also provide a powerful method of cluster 
detection. In this spirit, Schneider (1996) introduces the aperture mass measure as a way to systematically search for mass concentrations 
using weak lensing. In this approach, a weighted sum of image ellipticities is used as a proxy for the projected mass contained within an 
aperture. 
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Since gravitational lensing probes mass concentrations in a way that is independent of thieir baryonic content, shear-selected sam- 
ples can potentially provide a new and exciting view of large-scale structure. In fact, there have bee n several shear-selected cluster candi- 
dates that app ear to lack the characteristic galaxy over-den s ity and X-ray lumi nosity. (e.g.,[Erben et al. 2000; Umetsu & Futamase 20q3; 
iMiralles et"aUl2 002: Dah le et al.ll2003l : ISchirmer et al.ll2007l : iMaturi et alJlioOTh . The first detection of these so-called dark lenses was re- 
ported by Erben et abCOOOl) . Their initial analysis indicates significant tangential alignment roughly 7 arcminutes South of the cluster Abell 
1942. Assuming the presence of a mass concentration with roughly the same redshift as Abell 1942 {z ~ 0.2), they obtain a lower-bound 
mass estimate of M ~ 10^* M© inside of a sphere of radius r = 0.5/i~^ Mpc. However, a follow-up analysis bv lvon der Linden et al.l 
( I2OO6) using HST observations detects the tangential alignment with a lower significance, casting doubt on the hypothesis that the object is 
a true mass concentration. The HST data contains a larger number of distant galaxies and should therefore increase the significance of the 
detection in the case of a dark mass co ncentration. 

Among the most recent detections. lSchirmer et aljfeoOTh find 86 dark lenses out of 158 possible mass concentrations in a wide field with 
average galaxy density n ~ 12 arcmin"^. These objects, which show no detectable optical component, were identified using the aperture 
mass measure a nd a variant of it. Si milarly, using the linear filter of Maturi et al. ( 2005) to minimize spurious signals created by large-scale 
structure (LSS). lMaturi et al] ( |2007|) detect 7 dark lenses with no optical or X-ray component out of 14 identified lenses. The exact nature of 
these detections remains an open question. 

There are several possible explanations for the appearance of dark lenses in shear- selected surveys. One possibility is that these objects 
truly correspond to dark matter concentrations that are abnormally deficient in baryons. A significan t abundance of such obj ects is not 
expected and would require a rethinking of current structure formation scenarios jvon der Linden et al.ll200& : .Maturi et alj|2007l) . A second 
possibility is that dark lenses are spurious signals in a lensing map resulting from the alignment of intrinsic galaxy ellipticities or LSS 
projected along the line-of-sight. Both scenarios are likely to b e significant problems for weak-lensing surveys and have been investigated 
by se veral authors (Reblinskv & Bartelmann 1999; White et aL, 2002: Hamana et al. 2004; Hennawi & Spergel 2005; Pace et al. 2007; Fa^ 
l2007h . 

Another interesting possibility was proposed bv I Weinberg & Kamionkowskil ( |2002|) (WK2002). They suggest that dark lenses may be 
cluster progenitors that are not fully virialised. They argue that these objects should have a low galaxy over-density and X-ray luminosity 
compared to fully virialised clusters of the same mass. If some of these objects are sufficiently massive and over-dense to create a detectable 
weak-lensing signal, then they might have the same observable signatures as dark lenses. For clarity, we will refer to these objects as lensing 
protoclusters (LP). Using an analytic approach based on the Press-Schechter formalism, WK2002 estimate that 10 — 20 per cent of weak 
lenses should be dark LPs. 

Of course, it is entirely possible that the dark lens phenomenon is due to a combination of the above scenarios. Determining the extent 
to which each contributes to dark lens abundances, if at all, is an important step towards understanding the reliability of shear- selection 
techniques. In this paper, we explore the scenario proposed by WK2002 in further detail. We create a simple analytic model to investigate 
the likelihood that LPs can have the observational properties of dark lenses. 

The remainder of this paper is organized in the following manner. In section[2l we briefly review the aperture mass technique. In section 
[S] we present our analytic model of LPs. In section [4] we calculate the properties that LPs must have to be detected as dark lenses. We 
calculate the dark LP mass function in section[5]and explore what the abundance of dark LPs implies for weak-lensing surveys. Finally, we 
discuss our results and other potential dark lenses in section|6l 

In what follows, we assume a ACDM cosmology with parameters Q.m = 0.26, Ha = 0.74, Qb = 0.044 , H = WOh kms~^ Mpc~^ 
(with h — 0.72), n = 0.96, and as — 0.8, consistent with five-year WMAP results jDunklev et alj|2008l) . Unless otherwise noted, all 
distances and volumes are reported in physical units. 



2 THE APERTURE MASS MEASURE 

We use the aperture mass signal-to-noise ratio to define a detectable weak lensing signal in a shear-selected cluster survey. In this section, we 
briefly review the aperture mass technique (Schneider 1996). 

In order to realistically model a survey's sensitivity to mass concentrations, a functional form for the redshift distribution of background 

i ;alaxies in accordanc e with observations needs to be assumed. Following the analysis of WK2002, we assume a source redshift distribution 
Brainerdet"ai]| 19961) , 

with a mean redshift of 1.2 (zq ~ 0.8 and (3 = 1.5). In what follows we will find it convenient to isolate the source redshift dependence of 
all physical quantities by defining a redshift weight function (Seitz & Schneider, 1997. ; , Weinberg & Kamionkowsld..2002) . 

limz^^aoT,crit{zi;Zs) . Ecrit,oo(2i) 

Z{zs;zi) = H{zs-zi) = - -, (2) 

ljciit(zi; Zs) l^ciit(zi; Zs) 

where Zs and zi are the source and lens redshifts respectively. Edit is the critical surface mass density, and H is the Heaviside step func- 
tion. Note that the -dependence of the convergence, k — E/Ecrit, where E is the surface mass density, can be factored out by writing 
k{6, Zs, Zl) — 1^00 {9 , zi)Z{zs\ zi), where Kao{0,zi) = E(6')/Ecrit,oo(zO is interpreted to be the convergence for the case of reference 
sources at Zs — 00. Similarly, the shear can be written as 7(6, Zs, zi) = 700(6, zi)Z{zs\ zi). We will exclusively use the Zs -independent 
versions of the convergence and shear from here on. 
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The aperture mass provides a way to detect mass concentrations througii the distortion of lensed background images. It is defined as a 
weighted integral over the convergence, 

M^^iOo) = J d'e k{9 ) u{\e- e'o\), 0) 

where U{\9 — So |) is a compensated weight function centered at 6q, satisfying J dO U{9) = 0. As desired, Map can be expressed in terms 
of the tan gential component of th e shear relative to the direction 9o, 7t = — 5R[7(0 + do) exp{—2i<f))]. Using the so-called Kaiser-Squires 
mversion jKaiser&Squireslll993h . one can obtain 

Map(0o) = I d^e 'ytie- Oo) Q{\e^ eo\), (4) 



where Q{e) = ^ j^^ dx xU{x) - U{e). 

Equations l[3) and l|4j illustrate how the tangential shear field can be used as a measure of the integrated mass within an aperture. 
However, in practice, a discrete analog to equation Q in terms of a sum over individual image ellipticities is more useful: 

Map(e;) = i y Q{% - M) eu{e,; Bo). (5) 
n — ' 

i 

Here, n is the number density of galaxy images and eti{Oi \ 60) is the tangential component of the ellipticity relative to Oq. 
Setting 60—0 for simplicity, the expectation value of A/ap is given by 

(Map) = ^ ^Q(i0;i)7t.oo(e;), (6) 

i 

where {Z) = dzg Pz{zs)Z{za; zi). Note that in obtaining equation {6](, we have used {et) ~ (Z) 7t,oo in the weak lensing regime. In 
the case with no lensing, (Map) = 0, the dispersion of i\fap is obtained by squaring equation ^ and taking the expectation value. Assuming 
that tangential ellipticities of individual galaxies are uncorrelated, {en etj) = for i 7^ j, the dispersion is 

i 

where is the dispersion of intrinsic galaxy ellipticities. 

Averaging equations ^ and ^ over the probability distribution of galaxy positions and taking their ratio yields the ensemble averaged 
signal-to-noise ratio, 

S 2{Z)^ AO 6 {jt,^} (6) Q{e) 
N ^ Z ■ 



where Oout is the angular radius of the aperture and (74,00) (d) is the average tangential shear on a circle of angular radius 9. Using the 
Cauchy-Schwarz inequality the maximum signal-to-noise ratio is obtained by selecting Q{6) oc (7t) (6). This fact is intuitively clear; the 

signal is maximized by choosing the shear profile as its ow n weight functio n. 

In this paper, we use the weight function developed bv lSchirmer et alJ (l2004h (see also lSchirmer et al JzOOTI) . 

^ tanh(X/..) 

X/Xc 

where 

Here, X = O/Oout and Xc is a dimensionless width parameter. Smaller Xc results in more weight towards small radii. The above filter was 
developed as a computationally inexpensive alternative to using the NFW shear profile. It is designed to optimally detect NFW haloes in a 
wide field survey. Since our goal is to determine whether LPs are a significant contaminant in cluster surveys, the above filter will provide a 

more realistic estimate of their signal-to-noise. 

In what follows we use S/N = 4 as our shear-selection threshold (see lSchirmer et alll2007l for example). We adopt the fiducial values 
of n = 30 arcmiii^^ and cr,, — 0.3. In practice, a wide variety of aperture radii are used in order to detect clusters of different scales 
and at various redshifts. Since we aim to detect extended objects that are still in the process of collapsing, we adopt an aperture size of 
RovLt = Qo\itDA{zi) — 5 Mpc, where DAizi) is the angular diameter distance to the lens. This form is ideal for detecting objects with a 
scale of ~ 5 Mpc at any redshift. Note that 6'out becomes redshift dependent and corresponds to using 9out ~ 26, 14, 12, and 11 arcminutes 
for LPs at redshifts of 2: = 0.2, 0.5, 0.7, and 0.9 respectively. 
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3 ANALYTIC MODEL OF LENSING PROTOCLUSTERS 
3.1 The LP mass profile 

We qualitatively define LPs to be progenitors of cluster-scale fialoes tfiat are not fully virialised. Since haloes gain mass by accretion along 
their outskirts, it is natural to suspect that a protocluster consists of a small region in virial equilibrium surrounded by an infall envelope. 
We can therefore separate the protocluster into two distinct parts: the central virialised region (CVR) and the infall region (IF). In what 
follows, we define the virial radius rwi such that the average over-de nsity inside of ryir i s > 200 times the critical density of the universe at 
that epoch - a convention frequently used in N-body simulations. As lCuesta et al.l jlOOSl) point out, this definition does not provide a robust 
approximation to the true virial radius. We nonetheless adopt it for its simplicity and convenience in comparison to other works that have 
used this convention. Using different definitions of rvir does not change the main conclusions of this paper 
We model the mass density inside of the CVR with the ubiquitous Navarro-Frenk- White (NFW) profile, 

P{r)= . . , (11) 

where ps and are free parameters ( iNavarro et al.lll99^ll996l . n"997h . lTavio et al 1l l2008l) have shown that i ll It provides a good fit to density 
profiles within the virial radius. It is often convenient to characterize the above profile with the concentration parameter, Cvir = rv\v/rs. By 
integrating equation i ll It out to rvir and using = 200/9c(z) 47rr5i,/3, where niv is defined to be the virial mass and pc is the critical 
density of the universe, the concentration parameter can be related to ps through 

_ 200pc(^) 4^ 

3 ln(l + cvi.) - Cvir/(1 + CviO ■ 

The concentration d epends on the virial mass and redshift of the halo under consideration. In this paper, we use the form jSeliaklbOOd : 
iTakada & Jainll2002l) . 

where M, (z = 0) is the present day non-linear mass scale {5c{z = 0)/cr(A/,) = 1). As lTakada & Jairil l l2002h point out, th e halo model 
using the above form is known to be in better agreement with the non-linear matter power spectrum compared to other choices llSeliakl200d : 
ICoorav et al. 2000). 

We turn our attention to a quantitative description of the infall envelope surrounding the central regions of a LP. The outskirts are 
likely anisotropic due to the fact that accretion occurs along filamentary structures as observed in numerical simulations. Moreover, since 
the infall region consists of smaller haloes, we expect sub-structure to play an important role. However, the incorporation of both of these 
characteristics is beyond the scope of our analytic model. Instead, we content ourselves with developing a spherically symmetric profile 
describing regions beyond rvir. 

In what follows, we assume that the density profiles of LPs do not fall off as steeply as the NFW profile for r > rvir; we model the 
infall regions with a power law p ~ r~", where a < 3. This choice is motivated by several results. Using N-body simulations. 'Eke et al.| 
l ll998l) find a significant deviation from the NFW form at large radii, especially at higher redshifts (see Figure 10 of their pager). Secondly, 
the excursion set formalism can be used to show that the average density within infall regions falls off more slowly than r~^ lIBar kana 20040. 
Most recently, it h as been pointed out that the NFW form provides a poor fit outside of the virial radius jPrada et al.l2006l : ICuesta et al.„200a : 
iTavio et al.|[200^ . Using N-body simulations, iTavio et al.l j2008l) develop a density profile that better describes these regions on average. 
Their form closely approximates the NFW profile for small radii, but the instantaneous logarithmic slopes typically transition from ~ — 3 
at rvir to ~ —0.2 at lOrvir. They also find large variations in the density profiles of individual haloes beyond 2rvir. Hence, rather than 
determine a through dynamical arguments, we will explore what types of infall profiles are required to produce a detectable weak lensing 
signal in section|4] 

Following the above discussion, we model LPs with the density profile 

p{r) = { Jjrjrtr <r<R (14) 

R<r 

where rvir is the virial radius of the CVR and is a truncation radius, introduced to keep the mass of the profile finite. Throughout the rest 
of this paper, we will refer to equation l ll4b as the LP profile. 

Note that l ll4b is uniquely determined by four parameters: M (the total mass inside of the truncation radius R), rriv (the virial mass), 
zi (the redshift of the LP), and a (the logarithmic slope of the density profile outside rvir). Given m„ and zi, one can immediately obtain 
rvir using rriv = 200pc(z) 47rrvir/3. Equations l ll3b and l ll2b may then be used to obtain Cvir and ps respectively. Finally, po is obtained by 
imposing continuity at r = rvir. 



3.2 The surface mass density and tangential shear 

Let us define a lens-centered coordinate system, {n, r2, ra}, such that rs is along the line-of-sight to the lens center Our first task is to 
calculate the surface mass density, E(ri, r2) = drs p{r), of the LP profile in the 2 regions: the CVR and the outer infall region. 
Defining 
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Figure 1. Panel (a): the solid line shows the tangential shear of sources al Zs = oo induced by the LP profile. We assume parameters of A[ = 10^^ Mq, 
TUv = 5 X 10^"^ Mq, a = 0.5 and z; = 0.5. The dashed and dotted lines show the contributions from the CVR and infall region respectively. The dot-dashed 
line corresponds to a point mass of mass M. The tangential shear corresponding to the LP is the sum of contributions from the CVR and infall regions. Panel 
(b): the tangential shear due to the LP profile for values of o = 0.5 (solid), 1 (dashed), 1.5 (dotted), 2 (dot-dashed). All other parameters are the same as in 
(a). Protocluster density profiles that fall off more slowly with radius result in steeper up-tums in the tangential shear. 



{i:^-l)(l + Cvir) 
3(c2. -1)3/2 



(l_l2)3/2 



(l_j;2)3/2 



c2. -x2 



(3;2_l)3/2 



c2. -x2 



(l2_i)3/2 



a; < 1 



a; = 1 



x> 1, 



we obtain 

E(a;)LP = E(a;)cvR + S(a;)iF, 
where 

{2psrsf{x) X < Cvir 
a; > Cvir 

is the surface mass density of a NFW profile truncated at r^h and 



E(a;)iF = < 



y/cl,^ - x^ 2F1 



(15) 



(16) 



(17) 



I 



is the contribution from the infall profile, 



x < 

Cvir < a: < a;n 

X > XR. 



(18) 







Pi'r) — \ (r/r^ir) 

R<r 



r < rvir 
rvir <r<R 



(19) 



Here, 2F1 are hypergeometric functions, x = ( l/r^ ) -^/r J + projection of the coordinate vector in the plane perpendicular to the 

line of sight in units of r^, and xh = R/ra is the dimensionless truncation radius. 

The fact that equation il6\ is a sum of contributions from the CVR and infall regions will prove to be extremely helpful when 
quantifying the lensing contribution from unvirialised matter in section 1331 The -independent convergence is obtained from il6l using 
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Figure 2. Panel (a): the aperture mass signal-to-noise ratio as a function of aperture radius Sout ■ The solid, dashed, and dotted Hnes correspond to the LP profile 
and contributions from the CVR and infall regions respectively. We assume an aperture width parameter of Xc = 2, M = 10^^ Mq, = 5 X 10^'^ Mq, 
and a = 0.29 (chosen so the S/N = 4 for Bout = 14 arcminutes). Panel (6): same as (a) but with values of Xc = 0.5 (solid), 1 (dashed), 1.5 (dotted), and 
2 (dot-dashed). 



Koo = S/Scrit,oo- The Corresponding average shear profile, (7t,cx)) (0), can be calculated using (74,00) = Koo(^) — (d), where Roo{d) 
and (koo) (9) are the average value of the convergence inside and on a circle of angular radius 6 respectively. Owing to spherical symmetry, 
{koo) (0) = KooiO) and (74,00) {0) = 7t,oo(6') for the LP profile. 

We now investigate the tangential shear of sources at Zs = cxa induced by the LP profile. The solid curve in Figure[T](a) shows 7t,oo 
with lens parameters M = 10^^ M0, = 5 x 10^'^ Mq, zi = 0.5, and a = 0.5. The left and right vertical lines represent the virial and 
truncation radii respectively. Here, r^ir = 0.64 Mpc and R = 3.3 Mpc, corresponding lo 9 — 1.8 and 9.1 arcminutes respectively. Note 
that there are two kinks in the tangential shear. These kinks occur because il4\ is an idealized density profile, with sharp boundaries at r^ir 
and R. The dashed and dotted curves in panel (a) correspond to contributions to the shear from the CVR and infall region respectively. The 
LP tangential shear profile is the sum of these contributions. For reference, we also show the tangential shear induced by a point mass with 
M = 10^^ Mq (dot-dashed). The plot shows that the LP 

^t.oo is well approximated by the CVR 7t,oc inside of rvir- The LP 7t,oo rises 
outside of rvir, where the infall region contributes more to the shear. Outside of the truncation radius, the tangential shear is equivalent to the 
case of a point mass with the same mass M. 

In Figure[T](6), we vary the power-law index a for fixed values of M = 10^^ M©, m„ = 5 x 10^'^ M© and zi = 0.5. The solid, 
dashed, dotted, and dot-dashed curves correspond to a = 0.5, 1, 1.5, and 2.0 respectively. The rise in 7t,oo is most pronounced for flatter 
power laws. 



3.3 The aperture mass signal-to-noise ratio 

In this section we investigate the signal-to-noise properties of the LP. The solid line in Figure[2](a) shows the LP signal-to-noise ratio as a 
function of aperture radius for zi = 0.5, M = 10^^ M©, m„ = 5 x lO" Mq, and a = 0.26 (chosen so that S/N = 4 for 7?out = 5 Mpc, 
corresponding to 9out = 14 arcminutes). The left and right vertical lines correspond to the virial and truncation angular radii respectively. 
Since equation ([8} is linear in (74,00), the LP signal-to-noise ratio is the sum of contributions from the CVR and infall profiles. The dashed 
and dotted lines show the corresponding CVR and infall signal-to-noise ratios. Figure |2] illustrates that, for smaller aperture radii, the LP 
signal-to-noise is dominated by the CVR. On the other hand, the infall envelope makes a significant contribution in larger apertures. 

The left and right curves in Figure|2](6) show {S/N)cvr and {S/N)w using the same profile parameters as above. The solid, dashed, 
dotted and dot-dashed curves correspond to aperture width parameters of Xc = 0.5, 1.0, 1.5, and 2.0 respectively. Since lower values of Xc 
result in more weight towards smaller radii, the CVR contribution increases as Xc decreases. Conversely, the infall contribution is suppressed 
as Xc decreases. Note however that Xc has only a mild effec t on the LP sig nal-to-noise ratio. In practice, multiple values of Xc would be 
used to select clusters in a shear-selected survey (see Schirm er et al. Since apertures with larger Xc are more likely to introduce 

contamination from LPs, we choose Xc = 2 as our fiducial value from here on. 

Figure [3] shows the fractions (S'/iV)iF/(S'/iV)LP (solid) and (S'/iV)cvR/(S/iV)LP (dashed) as a function of m„ for M = 10^^ Mq 
and zi — 0.5. The power-law index a is varied to consistently satisfy (S/A^)lp = 4 for 9out = 14. The top axis shows the a required 
to satisfy this condition. We truncate both curves at q = 3 since the infall profile should not fall off more quickly than the NFW profile at 
large radii. As expected, the contribution from the infall region is greatest for smaller In this case, the density profile within the infall 
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Figure 3. Fractional contribution from tlie infall region (solid) and CVR (dashed) to the total LP signal-to-noise ratio as a function of virial mass . We 
assume a fixed total mass of M = 10^^ Mq and zi = 0.5. The power-law index a is varied so that the total LP signal-to-noise is 4 inside of an aperture with 
dowt = 14 arcminutes. The top axis shows the a satisfying this criterion for each rriv . Smaller virial masses result in a higher contribution from the infall 
region and flatter density profiles. 



envelope must fall off slowly with radius in order to meet the signal-to-noise threshold. On the other hand, larger virial masses result in a 
higher contribution from the CVR and steeper infall profiles. The plot shows that the infall region makes a significant contribution to the total 
signal-to-noise, even up to m„ ~ 1.5 x 10^* M0. 



4 LENSING PROTOCLUSTERS AS DARK LENSES 

In the last section, we described a simple analytic model for LPs that allows us to compute the shear and aperture mass signal-to-noise ratio. 
In what follows, we use the model to investigate the characteristics that a LP must possess in order to have the same observational signatures 
as a dark lens. 



4.1 X-ray luminosities and virial masses 

It is well known that a cluster's X-ray luminosity scales with its virial mass. Therefore, in order for a LP to be a plausible dark lens candidate, 
its virial mass must b e low enough to avoid detection in relatively deep X-ray searches. In this section, we use the semi-analytic calculation 
bv lNordetal.l (l200Sl to estimate the range of virial masses that a LP must have in order to be "dark". 

Given a sample's soft-band (0.1 — 2.4keV) flux threshold, iNordetaljJiooa) use the window function in equation (2) of their paper to 
model the fraction of virialised clusters detected (which we denote as /det) as a function of redshift. Figure |4] shows /dot for a flux limit of 
10~^* ergs cm~^. We assume a low flux limit here since a dark lens detection would likely be followed by a deep X-ray search. The 
solid, dashed, dotted, and dot-dashed lines correspond to virial masses of = 2, 4, 6, and 8 x 10^^ Mq respectively. Following lNord et"al] 
( l2008h . we assume a luminosity dispersion of oi = 0.59. At low redshifts {z < 0.05), nearly 100% of virialised haloes with the above masses 
are detected as X-ray sources. The percentage quickly declines with redshift. At z = 0.5, roughly 0, 1, 15, and 43 per cent of virialised 
haloes are detected with masses of 2, 4, 6, and 8 x lO'^^ Mq respectively. All detection fractions drop to nearly zero by z = 1. 

At a given redshift, we would ultimately like to estimate the maximum virial mass that a LP can have while still maintaining a small 
chance of being detected. For this task, we assume a fixed detection fraction /det and solve for the corresponding virial mass mjj as a function 
of redshift. Less than the fraction /det of virialised haloes are detected below this virial mass limit. The solid, dashed, and dotted lines in 
Figure |4](b) show the virial mass limit mjj for /det = 0.1, 0.3, and 0.5 respectively. As an example, less than 10 % of virialised haloes 
are detected below a mass of mjj — 5.5 x 10^^ Mq at z — 0.5. As expected, the plot shows that LP virial masses must be smaller at 
lower redshifts in order to maintain a significant chance of being undetected. From here on, LPs with a low probability of being detected via 
soft-band X-ray emission will be called "dark." 

Note that in the above cal culation, w e have assumed that soft-band X-ray luminosities of LPs roughly follow the mean scaling relation 
given by equation (1) of ,Nord et al.l feOOsh . Determining the extent to which this assumption is valid is difficult due to the effects of accretion 
and mergers on a cluster's X-ray luminosity. Using hydrodynamical simulations. iRowlev et alj J2004h find that the accretion of sub-clumps 
creates scatter in the L — niv relation by shifting clusters below the mean curve. They attribute this to the fact that while both the mass 
and luminosity increase as a sub-clump falls toward the center of the cluster, the temperature typically stays constant or decreases slightly. 
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Figure 4. Panel (a): the fraction of virialised haloes with a given mass detected via soft-band X-ray emission. We assume a flux limit of 10^^^ ergs cm~^ 
and luminosity dispersion of crj = 0.59. The solid, dashed, dotted, and dot-dashed curves correspond to halo masses of m„ = 2, 4, 6, and 8 X 10^^ Mq 
respectively. Panel (b): the virial mass limit m^J for a fixed detection fraction /dot = 0.1 (solid), 0.3 (dashed), and 0.5 (dotted). For a given redshift, less than 
/dct of haloes with mass < mjj are detected. 



Unfortunately, the task of adequately addressing the above issue is beyond the scope of our simple analytic approach. Further numerical 
studies are required to include these effects. 

It should also be noted that our discussion is restricted to cases where the mass growth rate of a LP is dominated by the accretion of 
smaller sub-clumps (ie that the CVR is the largest progenitor halo). In these cases where the sub-clumps are significantly cooler than the 
CVR, we expect the latter to be the dominant contributor to the integrated X-ray luminosity. This assumption may not be true for cases where 
the infall region contains a group that is of similar mass to the CVR. However, these objects should be morphologically different from our 
model, and therefore represent a different class. 



4.2 Aperture mass detection 

Armed with an appropriate range of dark LP virial masses, we turn our attention to shear selection. In this section we explore the physical 
characteristics of LPs that meet the aperture mass detection threshold of S/N — 4. 

In order to compare the densities of shear-selected LPs to virialised clusters, we calculate the minimum LP over-density required to 
meet the S/N threshold. The S/N of a LP is a function of M, m„, zi, and a. For a fixed M, m„, and zi, we solve for the power law index 
Q such that S/N = 4 inside of a lens-centered aperture. We then calculate the truncation radius R and the average over-density, Sf^, using 
1 + Sf^ = ZM / {A-KpE?'), where p is the mean matter density of the Universe. 

The solid, dashed, and dotted curves in Figure[5](a) show Sf'^ as a function of M at z; = 0.5 for m„ = lO", 5 x lO", and 
10^* M0 respectively. These masses correspond to /det ~ 0, 0.05, and 0.7 at z = 0.5. The dot-dashed line shows the virialisation threshold, 
^NL _ 200pc//O — 1 = 368, at z = 0.5. Panels (&) and (c) show the corresponding power-law index a and truncation radii R. We assume a 
signal-to-noise threshold of S/N — 4 and aperture radius of Rout ~ 5 Mpc for all curves. Figure[6]shows Sf^, a, and 7? as a function of z 
for M = 8 X 10^" (solid), 10^^ (dashed) and 3 x lO^'"^ (dotted) M©. Here, we assume a fixed virial mass of m„ = 5 x 10^^ Mq. 

In panel (a) of Figure[5]the minimum of Sf^ occurs when the truncation radius is similar to the aperture size. The aperture mass measure 
is most sensitive to overdensities with scales close to the aperture radius. The plots show that LPs must be increasingly over-dense to meet 
the S/N = 4 threshold as M decreases beyond the minimum, since the amount of mass within the aperture decreases. The over-density 
rises as M increases beyond the minimum because, as more mass is added to the regions outside of the aperture radius, the integrated shear 
inside of Bout decreases. In other words, although S/N increases at larger radii, the S/N within Oout actually decreases. Hence, a higher 
over-density is required to produce S/N = 4 as the scale of lens exceeds the aperture radius. 

Figure |5](fe) shows that LPs with a low probability of being detected via their X-Ray lu minosities must have infall regions with flat 
power-laws in order to meet the aperture mass detection threshold. We use the results of lTavio et a l. ( 2008) to qualitatively determine whether 
such objects are common in their high resolution N-body simulations. The fits to equation (12) of their pape r represent the density profiles 
obtained by averaging over all haloes in a given mass bin. For fixed parameters given in lTavio et alj ( |2008[) . we calculate the effective a 
required to reproduce the mass enclosed by their profile inside of lOrvir. Using a virial mass of 4.8 x 10^^ Mq and concentration 6.82, 
we calculate a value of q = 2.04, indicating that the flat profiles shown in Figure[5](&) are likely to be rare. In the next section, we will 
analytically estimate how rare these objects are using the excursion set formalism. 
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Figure 5. Panel (a): the average over-density required to produce an aperture mass signal-to-noise ratio of 4 inside of -Rout = 5 Mpc as a function of the total 
LP mass M. The solid, dashed, and dotted lines con'espond to my = 10^'', 5 X 10^^, and 10^'* Mq respectively. Also shown is the virialisation threshold 
at 2 = 0.5 (dot-dashed). Panels (b) and (c) show the infall power-law index a and truncation radius R required to satisfy the signal-to-noise condition. The 
aperture mass technique is most sensitive to LPs with R ~ Rout ■ 




Figure 6. Same as in Figure [5] shown as a function of redshift for a fixed virial mass of = 5 X 10^^ Mq. The solid, dashed, and dotted lines correspond 
to total masses of M = 8 X lO^'*, 10^^, and 3 X lO^^ Mq respectively. 



The plots also show that there are minimum and maximum detectable LP masses, which we denote as Mminirriv, z) and Mniax(»Tii,, z) 
respectively. For a fixed virial mass and redshift, Afmin(mi,, z) and iV/max(?ni,, z) are determined by where a — for constant S/N — 4. 
Using Q = as the minimum allowed power-law index ensures that we never consider cases where the LP density profile increases with 
radius. On the other hand, we could just as easily use some other non-zero minimum power-law index, amin- As Figure [5] (6) shows, the 
effect of choosing some Qmin > is simply to reduce the interval of detectable M. 

Finally, we note that it is the signal-to-noise requirement that determines a and R; we impose no dynamical constraints on these 
parameters. Therefore, not all of the values of a and R shown above may be physical. As an example, consider the dotted curve in Figure|5] 
(a), corresponding to = 10^* M©. It indicates that solutions exist for Af ~ 10^® M© that have extremely flat (q ~ 0) density profiles 
that extend out to i? ~ 6 Mpc. 

One way to eliminate some of the more extreme cases is to impose a higher value for amin- As Figures [5] (6) and (c) show, doing so 
both ensures that the density profiles fall off reasonably with r and removes the cases with the largest truncation radii. Ho wever, without any 
other dynamical arguments, choosing an appropriate Qmin is somewhat arbitrary. Given the large variations observed bv^ Tavio et al.l | |2008|) 
beyond the virial radius, extremely flat logarithmic slopes are possible, though the abundance of such infall regions has yet to be investigated. 
Since there is no obvious choice for Qmin, we will display results for multiple values in section [531 



5 ABUNDANCES 

5.1 The excursion set formalism 

The excursion set formalism ( Press & Schechte ill974l : lBond et al.ll99ll : [Lacev & ColJl99 3h was developed to infer the statistical properties 
of the non-linear density field using the framework of linear perturbation theory. Perhaps the most well known example is the derivation of 
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the halo mass function from linear theory, often referred to as the Press-Schecther (PS) mass function. However, one great advantage of the 
excursion set approach is that its basic framework can be applied to a wide variety of abundance calculations. Its versatility stems from the 
fact that, in principle, it may be applied to arbitrary linear over-density thresholds. These thresholds are analogous to the critical density 5c 
in the Press-Schechter halo abundance problem. In what follows, we will describe how the formalism may be used to obtain the abundance 
of a general object A defined by the linear over-density threshold 5a. In section we will fix 5a- 

Consider a point x in space within a realization of the matter density field at any early epoch in the Universe, before the growth 
of perturbations enters the non-linear regime. Rather than attempt to solve the non-linear evolution of density perturbations, we linearly 
extrapolate the initial density field to a later epoch using the growth factor from linear perturbation theory, D{z). The linear over-density 
at a point x at a later time z is simply given by 5{z, x) = 5{zi, yc)D{z) / D(zi), where Zi is the initial redshift. We then assume that the 
statistical properties of the true density field (for example, the halo abundance) at a given redshift can be inferred to a reasonable extent from 
the linearly extrapolated initial density field. 

For convenience, it is common practice to linearly extrapolate the initial density field to the present day. The over-density at a point x 
becomes redshift independent, but the over-density threshold becomes 5a{z) — 5a/D{z), where we have normalized D{z) to unity at the 
present day. In what follows we adopt this convention. Whenever it is necessary to discuss quantities that are not linearly extrapolated to the 
present day, we will note it in the text. 

The linear over-density around the point x is smoothed with a window function W{r; Rw) of scale Rw to obtain 



5{Rv 



d^rW{r-Rw)5{r). 



(20) 



Owing to the mathematical simplicity it affords, the most common choice for the window function is the fc-space top-hat window, defined by 



W{k;Rw) 



1 (k < R^) 
(k>R^). 



(21) 



One starts by smoothing the density field around x for large Rw (small k), and lowering Rw in increments. For each Rw, the variance 
of overdensities smoothed on this scale in an ensemble of density fields is calculated. 



S{Rw) = cj'^{Rw) = -r^ 



dk k^P{k), 



(22) 



where P{k) is the linear power spectrum. The set of points {S{Rw), 5{Rw)} traces out a trajectory parameterised by Rw in the {S, 5}- 
plane. In the limit that A.S O. lBond et al.l ( Il99lh showed that the probability density Q{S, 5) for a trajectory at {S, 5} satisfies 



dS 



2 ■ 



(23) 



We now turn to the case where there is an over-density threshold defining object A as discussed above. When 5{Rw) moves above or 
below 5a (depending on the particular application) at a scale S{Rw), the point x is assumed to be within an object of that scale. The goal 
then is to calculate the fraction of trajectories that cross 5a between the scales S and S + AS. Mathematically, this is realized by solving the 
diffusion equation J23b with absorbing barrier 5a. When the boundary condition Q(S, 5a) ~ and initial condition Q{So, 5) — 5o{5 — 5o), 
where So is the Dirac delta function, are applied to equation ([23}, Q{S, 5\So, 5o) d5 represents the probability that a trajectory starting at 
{So, 5o} obtains an over-density between 5 and 5 + d5 at S without having crossed 5a. The fraction of trajectories that cross the threshold 
at or prior to S{Rw) is given by the complement of Q, 



F{S,5a\So,5o) = 1 



Q{S,5\So,5o) d5. 



(24) 



Equation ( 124b . which applies only to the case where up-crossings are of interest, represents the fraction of mass within objects A with mass 
greater than M[Rw). Hence, the differential fraction of mass within objects A is given by 



fs(S,5A\So,5Q) = 



dF{S,5A\So,5o) 
AS 



85 



(25) 



where the last equality was obtained by using equation ( 123b . Equation l l25b is often referred to as the first-crossing distribution. Taking 
{5*0, 5q} — {0, 0} (ie that the density field approaches the mean when smoothed on large scales), the mass function of object A may be 
obtained from 



n(M) dM = 



M 



dF{S,5A) 



dS 



dS 
dM 



dM. 



(26) 



Note that up until this point, we have assumed the use of the fc-space top-hat window function. However, the main disadvantage of 
equation ( 121b is that both the volume and mass within W{k; Rw) is not well defined. To overcome these problems, it is common practice to 
derive equation ( 126b using the fc-space filter, but at the end replace the variance ^(fc) by the real space top-hat relation, 



S{Rv 



27r2 



dk k^P{k) 



3 sin(fc7?i^) — ZkRw cos(fc_R„) 
{kR^^ 



(27) 
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Figure 7. The linear over-density of a LP required to create S/N = 4 inside of -Rout = 5 Mpc versus the variance of the linearly extrapolated density field, 
t7^(Af). From top to bottom, the solid curves show the exact linear over-density for m„ = 10^^, 5 X 10^"^, and lO^** Mq respectively. The dashed lines 
show the corresponding piecewise approximations (see appendix|A). Note that the ordinate is not linearly extrapolated to the present day. 



In this case, the mass within the Lagrangian radius Rw is given hy M — p 47r_R^/3. Therefore, in the final relations that we derive, we are 
to interpret the quantity as the variance of the smoothed linear over-density inside of the Lagrangian radius Rw - Similarly, the smoothed 
over-density 5 is to be interpreted as the average over-density inside of Rw- 



5.2 The mass function of dark LPs 

In this section, we aim to estimate the abundance of dark LPs using the excursion set formalism Unfortunately, the formalism cannot be 
used to calculate the abundance of objects with density profile il4t . Instead, we content ourselves with the more modest goal of considering 
overdensities that are large enough to create a weak lensing signal, but are unlikely to be detected via their X-Ray luminosities. Our approach 
contains two steps: Step A) we add up the fraction of mass within overdensities above a lensing threshold. For this we use the overdensity 
barriers obtained in section 1431 Step B) we multiply by the fraction of these overdensities that have a low probability of being detected by 
their X-Ray luminosities. 

We begin our discussion with Step A. Since the excursion set formalism is based on linear perturbation theory, our first task is to convert 
the Sf^ barriers of section [42] to linear overdensities, 5i, using the spherical collapse model. Figure [7] shows the results of this conversion 
as a function of the scale S = cr'^(Af). M is shown on the top axis (Recall that 5 is a monotonically decreasing function of A/). Here we 
assume zi — 0.5 and Qmin = 0. From top to bottom, the solid curves correspond to the exact 5i for ttIu — 10^"', 5 x 10^'', and 10^* 
respectively. The dashed curves show the corresponding piecewise approximations discussed below and in appendix lAl 

Note that each barrier shown in Figure [7] is scale-dependent and cannot be expressed as an analytic function of S- Obtaining exact 
analytic solutions for their first crossing distributions is therefore impossible. For a fixed m„, we address this issue by approximating the 
barrier as two lines, with the form of equation dAlb . In appendixlAl we show that a solution for the first crossing distribution with an absorbing 
barrier of this form can be reduced to quadrature. In what follows, we will use the approximate first crossing distribution fs{S, Si), given by 
equation lA4l 

We make two additional notes about the over-density barriers in Figure]?] First, 5i is a relatively weak function of rriv, particularly for 
smaller values of m^- In fact, below mv — 10^'^ Mq the amplitude of 5i is virtually independent of m^. Secondly, for two virial masses m" 
and mj^, where > m', any trajectory that crosses 5i{m\,) must also cross Si{m^) at a larger total mass scale. We therefore assume that 
fs [S, 5i{mv)] yields the approximate fraction of mass in lensing overdensities that potentially correspond to dark LPs with virial mass less 
than rriv- However, many of the lensing overdensities satisfying 5 > 5i contain large sub-haloes. Following section l4n if the sub-halo masses 
are large enough, then they have a high probability of being detected via their X-Ray luminosities. Hence, their host overdensities would not 
satisfy the "dark" criterion. Our goal, then, is to calculate the fraction of overdensities whose sub-haloes do not exceed a maximum probability 
/dot of being detected via X-Ray luminosity (Step B). Using the results of section 1431 we can map this probability to a maximum allowed 
sub-halo mass, mjj. Put in another way, our goal is to obtain the fraction of mass contained in trajectories without "nearby" trajectories that 
cross the virialisation threshold at mass scales greater than mjj. 

In what follows, we adopt the notation Sv = cr^(m„). We also denote the linear virialisation threshold as Sy, w hich is obtained by 
applying the spherical collapse model to 5^^ = 200 p c {z) / p{z) — 1. Using N-Body simulations iCasas-Miranda et alj JloOZ) find that the 
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0.5. The solid and dashed Hnes correspond to /dot = 0.1 (m^ = 5.5 X 10^ 
cm^^ is assumed. The dotted line corresponds to /(j^t 
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Figure 8. Panel (a): the dark LP mass function at zi 

(8.5 X 10^^ M©) respectively. A soft-band X-Ray flux limit of 10"^** ergs s"^ ^^,.-2 
limit of 5 X 10~l'* ergs cm~^ (m^ = 1.4 X 10^* Mq). We assume a, 

function at the same redshift. Panel (6): the probability that a 5; over-density has zero sub-haloes with mass between mjj and M, shown for the same mjj 
values above. Since this probability is low in the mass range under consideration, the abundance of dark LPs is suppressed. Panel (c) : the infall power-law 
index required to satisfy S/N = 4 for = mjj. This panel serves as a reference for the types of LP profiles required to create a weak lensing signal in the 
mass range shown. 



probability, Pv(Nh, m\M, 5), of having Nh sub-haloes with mass greater than m within an over-density 5 of mass M is well described by 
a Gaussian with mean number 

I'M 

{N){m\M,S)= / dm^. N{mv,Sv\M,5) (28) 



-mi 



and variance 

VariNh)^ {l + AD^{z)S) I I dmidm2N{mi,3,\M,5) N{m2,S,\M - mi,5') + {N) - {Nf , (29) 

where 

, {S^ - 5) 

=dv - -, ijjT (30) 

1 — (mi/M) 

and 

N{m^,S^\M,S)dm^ = —fs^{Sv,S^\S,S)dm^ (31) 
ami, Triv 

is the average number of virialised sub-haloes with mass between m„ and m„ + dm„. The second term in the prefactor in equation ( 1291 ) 
is a phenomenological term accounting for clustering effects. For the constant A, we use 0.05, which was calibrated to simulations by 
ICasa s-Miranda et al.l J2003) . 

The probability that a 5i over-density has zero sub-haloes in the mass range (mjj, M) IS given by Pv(0,m^|A/, 5i). Therefore, the 
fraction of mass within dark LPs whose sub-haloes have a probability < /dot of displaying detectable X-ray emission is fs \_S, 5i{m'^)^ x 
Pv{0, m^\A4, Si), and the mass function is given by 

fs[S,5i{mt)] X Pv{0,mt\M,3i) (32) 



nD{M,z) = 



dM 



The solid and dashed lines in Figure[8] (a) show the dark LP mass function at — 0.5for/dct ~ 0.1 and 0.5, corresponding to sub-halo 
mass limits of — 5.5 and 8.5 x 10^'' Mq respectively. We assume the fiducial soft-band flux limit of 10~^* ergs cm~^. The dotted 
line shows the mass function for /dot ~ 0.1 and a flux limit of 5 x 10^^ ergs s~^ cm' minimum allowed 

power law index of Qmin = is assumed for all curves. The only notable effect of changing Qmin is to change the domain of detectable dark 
LP masses. As figures[5](a) and (b) show, increasing amin truncates the low and high mass ends of detectable overdensities. For reference, 
we show the lSheth & TbrmeiJ j 19991) (ST) halo mass function (dot-dashed) for z = 0.5. 

Panel (6) shows Pv [O, m'^\M, (5;(m^)] for the same values of m'^. The mass function of dark LPs is suppressed because the probability 
of finding a 5i over-density without sub-haloes large enough for X-ray detection is small. In addition, the large-mass end is suppressed in 
two ways: 1) The shape of the over-density barriers in Figure|7]make it less likely for trajectories to cross at higher mass scales. 2) Panel 
(6) shows that as the total mass M increases, it becomes less likely that the over-density will contain zero sub-haloes in the mass interval 
(m^, M). Panels (a) and (&) also show that increasing mjj through either /dot or the flux limit yields only a mild increase in the amplitude 
of the mass function. Up to ~ A//2, the probability factor Pv [O, m'^\AI, Si (mjj)] remains low for increasing due to the decreasing 
variance ( |29b . 
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Figure 9. Panel (a): the differential number count of dark LPs for a hypothetical weak-lensing survey with background galaxy density n = 30 arcmin"^ 
and intrinsic ellipticity dispersion cTf = 0.3. We assume a fixed /^ct = 0.1. The solid, dashed, and dotted curves correspond to a^\^ = 0, 0.2, and 0.4. 
The dot-dashed curve shows the estimated differential number count of clusters. Panels (b) shows the same but with fixed Qniin = and /jgt = 0.3 (solid), 
0.5 (dashed), and 0.7 (dotted). Panel (c): the fraction of lenses that are dark LPs, Fp, as a function of redshift. The dashed curve shows Frt using our 
redshift-dependent cluster mass threshold M^i{z). From bottom to top, the solid lines con'espond to constant M^i of 8 X 10^"^, 1,2, A, 6, and 8 X lO^** Mq. 



Panel (c) of Figure[8]shows the power law index a required to create S/N = 4 for the LP profile. The purpose of panel (c) is to provide 
a reference for the types of density profiles required to create a detectable weak lensing signal in this mass range. Note, however, that there 
is no direct relationship between a and the overdensities counted using the excursion set procedure in this section. This is a limitation of 
analytic approach taken here. The excursion set formalism does not yield information on the mass profiles of 5i overdensities. It is therefore 
impossible to rigorously quantify the aperture mass S/N within the formalism. However, we argue that by selecting objects with an adequate 
over-density to create a weak lensing signal, and correct sub-halo structure to avoid X-ray detection, we obtain a reasonable estimation of 
dark LP abundances. 

Finally, we note that we have not taken into account all of the trajectories that may correspond to dark LPs. Some trajectories may 
obtain 5 > 5; at a scale S = Si (see appendix |A]i, corresponding to the maximum detectable LP mass Mmax- A fraction Pv of these 
trajectories correspond to dark LPs with mass Mmax. However, since we have shown that large-mass LPs are extremely rare, we can neglect 
these trajectories with little consequence. 



5.3 Dark LP counts and weak lensing surveys 



Here we investigate the abundance of dark LPs and its consequences for future shear-selected cluster surveys. For a fixed /dot the number 
counts of detectable dark LPs per unit steradian, per unit redshift interval, is given by 

„ax(z) 

nDiM,z)dM (33) 

i„(^) 



dzdQ dzdO, 

where no is given by equation i32l . Here, 

dV _ c {l + zfDA{zf 
dzdfl ^ Ho ^Q^{1 + zY' + nA 

is the comoving volume element, where c is the speed of light and Ho is the present-day Hubble parameter. 
Similarly, the differential cluster number count is 



(34) 



dNci{z) dV 



dzdO, dzdfl 



nh(M,z)dA/, (35) 



where, Uh is the halo mass function, and Aid is the redshift-dependent mass detection threshold of the survey. We use the Sheth-Tormen 
mass function for Wh. To estimate Md for a hypothetical survey using the filter ^ with n — 30 arcmin^"^ and — 0.3, we use a NFW 
profile truncated at the virial radius. At a given redshift, we solve equation ^ for the mass that yields S/N — 4. Since a wide range of dout 
and Xc would be used in practice, we set Rout equal to the virial radius and Xc equal to the cluster concentration in order to estimate the 
lowest detectable mass. Using this method, we obtain detection thresholds of M^i = 8.2 x lO", 1.1 x lO", 1.4 x 10^", and 1.9 x 10^* Mq 
ai z = 0.1, 0.2, 0.3 and 0.4 respectively. 

The solid, dashed, and dotted lines in Figure|9](a) show the differential dark LP number count for Qmin ~ 0, 0.2, and 0.4 respectively. 
We assume a fixed /dot =0.1. This means that all LPs under consideration have sub-haloes with less than a 10 per cent chance of being de- 
tected via soft-band X-ray emission, assuming a flux limit of 10~^* ergs cm^^. The plot shows that most dark LP detections correspond 
to objects between z ~ 0.1 — 0.6. The dot-dashed line shows the differential cluster number count obtained using the AIci{z) described 
above. 
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Note that the most efficient lenses are located at 2: ~ 0.5 (roughly halfway between the observer and the peak of the source distribution). 
However, geometry is not the only factor affecting the differential number counts. For both dark LPs and clusters, the growth of structure leads 
to an increase in the comoving number density at lower redshifts. On the other hand, the differential number counts are suppressed in this 
regime due to the comoving-volume element. At high redshifts, the lack of LSS and higher detection limits are responsible for suppressing 
the number counts in both cases. The dark LP number counts decrease as Qmin increases because the low and high mass LPs with the flattest 
power-laws are cut out (see Figure [5]l. 

The solid, dashed, and dotted curves in panel (6) correspond to /dot = 0.3, 0.5 and 0.7 respectively. We assume a fixed amin = 0. 
Raising /dot increases the dark LP number counts since overdensities with larger sub-haloes are included. However, panel (6) shows that this 
is a mild effect; the results are relatively insensitive to the choice of /dct. 

Panel (c) shows the fraction of lenses that correspond to dark LPs, 

/^'■""''ni,(M,2) AM 

Fd{z) = (36) 

/ 1^2^ no (M, z) dM + /^^^_ nh(M, z) dM 

for fixed Qmin ~ and /dct = 0.5. From bottom to top, the solid curves correspond to constant M^i = 8 x 10^'^, 1, 2, 4, 6, and 
8 X 10" Mq. The dashed curve was obtained using the redshift-dependent AIci{z) described above. Note that in the denominator of i36l 
we have neglected the shear selected counts due to other potential causes of dark lenses (ie LSS and allignment of intrinsic ellipticities). 
Including these detections would decrease Fd- We will discuss spurious signals in section |6] Panel (c) shows that even in our worst case 
estimate, where sub-haloes in LPs have up to a 50 per cent chance of displaying observable X-ray emission and flat infall profiles are allowed, 
dark LPs only make up < 0.4 per cent of lenses at at any given redshift (see the dashed curve in Figure|9](c)). 

By definition, the redshifts of dark LPs are unobservable. Hence, weak-lensing surveys will only be sensitive to their cumulative number. 
We calculate the all sky number of dark LPs and clusters by integrating equations ([33} and J35b for Ofmin ~ and /dct =0.5. We obtain 
total dark LP and cluster numbers of 410 and 108, 303 respectively. We estimate that the former make up < 0.4 per cent of the total number 
of lenses. It is therefore unlikely that they will be r esponsible for a significant fraction of dark lens detections in future surveys. Although 
we restrict our analysis to the aperture mass filter o f ISchirmer etal] j2004 [2007). this result holds for other filters as well. A filter that adds 
more weight to its outer regions would decrease the over-density required for a LP to create a weak lensing signal. However, the main factor 
suppressing dark LP abundances - the probability of finding such an over-density without large sub-halos - remains small. We emphasize that 
the dark LP abundance is generally small because it is highly unlikely to find an over-density large enough to create an observable lensing 
signal with sub-haloes small enough to escape X-ray detection. 

In accordance with current shear-selected samples, we perform ca lculations using n = 24 arcmin"^ and = 0.48. These parameters 
are more representative of the deepest exposures in surveys to date (see lSchirmer et al.l2007l for example). In this case we find that dark LPs 
cannot be detected by their weak-lensing signal. Put in another way, there are no protoclusters with small enough virial masses that meet the 
S/N = 4 threshold. In order to create a detectable weak-lensing signal, the virial masses have to be larger. However, this means that the 
probability of X-ray detection increases. Hence our model indicates that current shear-selected samples are not likely to contain a significant 
number of LPs that will escape deep ( > 10^^^ ergs s^^ cm^^) soft-band X-Ray searches. 



6 DISCUSSION 

We have developed an analytic model to determine whether LPs can account for a significant fraction of dark lenses. In our model, a 
protocluster consists of a small virialised central region surrounded by infalling matter. A dark LP corresponds to a cluster-scale mass 
concentration with group-sized W^^ M©) virial mass. The small virial mass results in a low probability of X-ray detection, while 
the large total mass yields a high aperture mass S/N. As initially suggested by WK2002, these objects can potentially share the same 
observational properties as dark lenses. 

For the LP mass distribution we used an idealized model consisting of an NFW profile inside of the virial radius and power-law, 
p ~ extending from the virial radius to the truncation radius. In this case, the total S/N is the sum of contributions from inside and 
outside of the virial radius. In the case of small virial mass, the S/N is dominated by the infall region. 

Dark objects in shear-selected samples would likely be follow ed up with deep X -Ray searches. In order to quantify the likelihood for 
a LP to display detectable X-ray emission, we used the analysis of iNord et al.l ( l2008i) . We found that LPs with a low probability of being 
detected via their X-Ray luminosities (or equivalently low virial masses) must have large total masses (M ^ IO'^^Mq) an d a-values ^ 1 to 
meet the aperture mass detection threshold. Such infall regions may exist given the recent findings of Tavio et al] ( |2008|) . who showed that 



the density profiles of haloes deviate from the NFW form beyond the virial radii, and display considerable scatter. The abundance of these 
objects in N-body simulations has yet to be investigated. Objects with the above characteristics would display rising (74) at larger radii. A 
comparison of our results with the shear profiles of detected dark clumps is difficult due to the fact that their redshifts are unmeasurable by 
definition. Hence, it is impossible to determine whether features in the shear profile occur at the appropriate radii. 

We have used the excursion set formalism to calculate the abundance of dark LPs. In our approach, the number density of mass 
concentrations that are sufficiently overdense to meet the S/N threshold is multiplied by the fraction that are unlikely to be detected in the 
soft X-ray band. This subset of objects contains zero virialised sub-haloes with a probability > /det of being detected via X-Ray emission. In 
most cases of interest this fraction i s extremely small, resulting in a suppression of dark LP abundances. These results appear to be consistent 
with the average profiles derived in lTavio et al.l ( 1200 Sl) . which indicate that infall regions typically do not contain enough mass to create dark 
lenses. 

In section l53] we compared the differential number counts of dark LPs to ordinary clusters in a hypothetical shear-selected survey with 
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source density n = 30 arcmin ^ and intrinsic ellipticity dispersion CTe = 0.3. In both cases, we found thiat most detections originate from 
objects at 2; ~ 0.1 — 0.6. The dark LP number counts are generally 2 — 3 orders of magnitude smaller than the cluster number counts. 

We varied the minimum allowed power-law index, amin, to simulate scenarios in which flat infall power-laws are dynamically unlikely. 
We found that dark LP abundances are highly sensitive to ctmin, dropping rapidly with increasing amin. If infall regions typically fall off 
steeper than r^^, then lensing contributions from the outskirts of dark LPs may be insufficient to meet the detection threshold. We also 
varied /dot to explore the remote possibility that objects with large virialised sub-haloes could be detected as dark lenses. We found that the 
differential number counts are relatively insensitive to /dct, varying only by ~ 20 — 30 per cent between /dot = 0.3 and 0.7. 

Finally, we have calculated the fraction of lenses that correspond to dark LPs. We found that they constitute < 0.4 per cent of lenses at 
any given redshift. Moreover, dark LPs account for ^0.4 per cent of the total number of lenses in our hypothetical shear-selected survey. 
We therefore concluded that dark LPs are too rare to be considered a plausible dark lens candidate. 

Our approach adds to initial work by WK2002 in several important ways. The first is our use of the density profile ( I14t . which provides 
a physical model for lensing protoclusters that takes into account devi ations from th e NFW fo rm beyond the viri al radius. Such deviations 
have been recently pointed out in high resolution N-body simulations bv lCuesta et alj|2008l) and lTavio et alJjToOSh . An additional advantage 
of iI4b is that we are able to quantify the S/N contributions from virialised and unvirialised matter. In contrast, the model of WK2002 does 
not allow one to quantify the lensing contribution from the central regions that meet the virialization over-density threshold. Hence, in their 
approach it is possible to consider cases where the S/N is dominated by the virial mass of the LP. These cases typically occur when the total 
over-density of the LP is close to the virialization threshold. Note that this difference is one of the reasons that our approach generally yields 
higher LP masses compared to the WK2002 results. By forcing the virial region to be smaller in order to simultaneously minimize its lensing 
contribution and avoid X-ray detection, larger total masses are required to meet the S/N threshold. 

We also point out that the fraction of mass in lensing overdensities does not correspond to the fraction contained in dark objects. Many 
of these overdensities contain large virialised sub-haloes. In practice, these cases would corr espond to true cluster detect ions since follow up 
X-ray searches would be sensitive to these sub-haloes. By incorporating the halo statistics of lCasas-Miranda et al.l ( l2002h . our approach only 
counts overdensities without large virialised sub-haloes. This key difference accounts for the lower dark LP abundances that we obtained 
compared to WK2002. 

Following t he wo rk of iReblinskv & BartelmamJ ( Il999h : IWhite et al.l ( l2002h : iHamana et al] ( l2004h : iHennawi & Spergell ( l2005h : 
IPace et alj ( l2007h : IpaJ J2007h , it is more likely that dark detections will correspond to false peaks resulting from: 1) LSS along the line- 
of-sight. In this case, the S/N is due to projected mass; it cannot be associated with a single isolated structure. 2) the random or correlated 
alignment of intrinsic galaxy ellipticities. These alignments alone can lead to spurious detections, especially in shallow surv eys. However, 
it is also possible they can boost peaks that correspond to smaller mass concentrations dvon der Linden et alj|2006l ; lFai]|2007h . Owing to an 
artificially high S/N, these detections can be misinterpreted as dark lenses. 

Using ray - tracin g through stacked snapshots of cosmological N-Body simulations, IPace et alj j2007h tested the performance of the 
ISchirmer et alj j2004h filter used above (referred to as O APT in their paper). By removing individual lens planes of haloes that may be 
associated with a particular S/N peak. [Pace et al.l ( 120071) were able to separate true detections from spurious ones. A true detection is 
associated with a cataloged cluster in the N-body simulation; a spurious detection remains when lens planes of individual candidates are 
removed. They found that the OAPT filter yields spurious detection fractions < 20 (25) per cent at S/N — 4 for source redshifts of 
Zs = 1 (2). For larger aperture sizes, this fraction decreases only mildly at higher S/N thresholds. In addition, they point out that these 
spurious detections are indistinguishable from true detections in a S/N map. Therefore, it is likely that LSS accounts for at least some of the 
dark lenses reported in the literature. 

Note that a dark LP would likely be counted as a "true" detection with the algorithm of lPace et al.l ( l2007h . The removal of the lens plane 
containing the dark LP would significantly diminish the signal observed in the S/N map. In addition, since the CVR would be cataloged 
in the N-body simulation, the detection might be associated with this small -mass halo . There fore, it would be instructive to determine what 
causes the lensing enhancement of small-mass detections in studies such as lPace et alj j2007h . 

Finally, we point out that the intrinsic galaxy ellipticities were randomly oriented in Pace et al. (2007). The effect of correlated alignment 
of intrinsic ellipticites on the number of false detections was not taken into account. As iFan, C2007,) points out, galaxy formation is sensitive 
to the local environment. One would therefore expect the orientation of a galaxy to at least be correlat ed with its closest neighbors. Such 
alignments can increase the number of false detections in convergence K-maps significantly, "paiil ( |2007|) showed that including this source 
of noise can increase the likelihood of false detection due to intrinsic ellipticities in a given field. This increase can affect whether intrinsic 
ellipticities can be ruled out in a dark lens detection. Future numerical studies on false peaks in weak lensing surveys should investigate this 
important possibility. 

While the analytic model presented in this paper provides important insight into why dark LPs should be extremely rare, it is limited 
by several key issues. The first is the simplistic density profile JI4b . which neglects the effects of anisotropy and substructure on the S/N. 
A more detailed analysis should incorporate these properties, which are expected to have a significant effect on the lensing signal. Secondly, 
since the excursion set formalism does not yield any information about the density profiles of individual trajectories, it is impossible to 
rigorously determine whether objects meet the S/N threshold. The best we can do in our analytic approach is to assume that objects above 
the derived over-density threshold can be detected. In addition, we have used a simple model for the X-ray luminosities of protoclusters in 
order to estimate the probability of detection. In reality this is a highly complicated problem with many caveats that can only be addressed 
numerically. Lastly, our model does not include galaxy overdensities. We assume that objects with low virial masses also escape optical 
selection. Future studies should focus on whether LPs display low galaxy overdensities as well. Each of the above issues would be ideally 
addressed in a high-resolution N-body simulation containing a baryonic component. Our model provides a starting point for more detailed 
investigations on the characteristics of simulated protoclusters and their impact on shear-selected samples. 
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APPENDIX A: THE FIRST CROSSING DISTRIBUTION OF THE WEAK LENSING BARRIER 

In this section, we obtain the probability of piercing the absorbing barrier, 
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BaS + Aa 

BtS + {Ba - Bb)S2 + Aa 



Sl<S<S2 
S2<S< 53, 



(Al) 



at a scale 5. The quantity Si represents an approximation to the minimum hnear over-density required to create an aperture mass signal-to- 
noise ratio of 4. Here, 5i = cr'^(A'/nmx) and 53 — cr^iMmin) are the mass scales corresponding to the maximum and minimum detectable 
LP masses discussed in the section|4] 52 is the mass scale at which 5i transitions from slope Ba to Bt- 

Fortunately, the problem can be greatly simplified by the fact that a trajectory beyond a given scale 5 is independent of the path leading 
up to 5. We do not want to count trajectories that simultaneously cross the virialization threshold, 5v, at 5 < 5i and 5; at 5i < 5 < 53. 
These correspond to lensing overdensities within larger collapsed objects. To avoid this problem, we use 

iL 

25i 



Qps{Si,Si) 



exp 



exp 



25i 



(A2) 



as the probability density for starting at the origin and ending at { 5i , (5i } . Hence, for example, the probability density for a trajectory starting 
at the origin with intermediate and end points of {5i,(5i} and {5,5} respectively, where 5i < 5 < 52, is (3ps(5i, (5i)(5"(5, (5|5i, 
Here, (5, (5|5i, 5i) is the conditional probability density for trajectories starting at {5i, where 5i < 5 < 52 • The superscript denotes 
that the barrier parameters in the regime Si < S < S2 are to be used. The general form for this probability density is given by equation l |B5l l 
in appendix IBI Similarly, if 52 < 5 < 53 then the probability density is Qps{Si, 5i)Q°(52, d2\Si, 5i)Q''{S, 5\S2,32), with the additional 
intermediate point {52, 52}- Summing over intermediate points yields 



Q{S,S) = 



dSi Qps{Si,Si)Q''{S,S\Si,Si) 



5i < 5 < 52 



(A3) 



BaS2+Aa rBaSi+A 



I- 



dS2 dSi Qps(5i,5i)Q"(52,<52|5i,5i)Q''(5,5|52,52) 52 < 5 < 53 



The first crossing distribution fs{S, Si) of dAlt for trajectories starting at the origin is 



fs{S,Si) = 



dS 

BaSl+Aa 



dSi Qps{Si,Si)f^{S,5i\Si,Si) 



(A4) 



5i < 5 < 52 



X 



BaS2+Aa fBaSi+A 



d<52 d5i Ops (5i , <5i ) 0" (52 , <52 1 5i , 5i ) /|(5, 5, 1 52 , 82 



52 < 5 < 53 



where /s(5, 3i\Si,5i) and /I (5, Si\ 82,52) are obtained from equation JB7t . 

To illustrate the characteristics of l lA4t . we compare it to the PS first crossing distribution at z = (depicted with crosses) in Figure 
lAll (g). For the latter we use 5„ — 1.63, which is obtained by applying the spherical collapse model to the virialisation threshold 5^^ — 
200pc/p — 1. The solid, dashed, and dotted curves correspond to iA4\ with barrier parameters Ba = —Bb = 0.001, —0.1, and —0.5. For 
all curves, we assume 5i = 0.1, 52 ~ 1.0, 53 — 5.0 and Aa = Sv Panel (6) shows the corresponding barriers. 

The kinks in panel (a) correspond to the mass scale 52 at which the two linear barriers meet in iAH . The solid curve shows that, in the 
limit where iAl\ is approximately constant with a value of ~ 5^, equation iA4l is equivalent to the PS form. As the absorbing barrier dips 
down, it is more likely for trajectories to be absorbed in the 5i < 5 < 52 regime. In this case, the number of available trajectories to pierce 
the S > S2 side of the barrier is depleted. This effect can be observed in panel (a) as a decrease in the first-crossing probability for S > 82- 



APPENDIX B: CONDITIONAL FIRST CROSSING DISTRIBUTION FOR THE LINEAR ABSORBING BARRIER 

In this section we obtain the first-crossing distribution of the linear absorbing barrier S = BS + A for trajec t ories starting at {5o,(5o}- 
Solutions to the problem in which trajectories start at the origin can be found in IShetlj (Il998l) andlMcOuinn et al.l |2003) . 

For completeness, we summarize the diffusion equation approach taken by iMcOuinn et al j2005l) to obtain the general solution for 
Q(5, S). The probability density obeys the diffusion equation 

dQ _ 1 d'^Q 
dS ~ 2 dS^ 

with boundary condition Q — Ofoi 5 — BS + A. We utilize the linear transformation y — B{5 — BS) and a; = 5 — 5o to obtain 

dQ^B^t^ + B'^ (B2) 
dx 2 dy^ dy 

with Q{x,y — BA) = 0. Assuming a solution of the form Q{x, y) = f{y)g{x), the problem is reduced to solving two ordinary differential 
equations, g' = \g and {B'^ /2)f" + B^f' = A/. The general solution can be written as an integral over the parameter A, 



(Bl) 



Q{x,y) 



-B'^/2 



+ 



-B2/2 
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d\ h{X) exp'''+^" (ci(A) exp'"^ +C2(A) exp"*"^) 
dX h{\) exp-^+^" (ci (A) exp"'" +C2(A) exp"'"") , 
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Figure Al. Panel (a): comparison of equation jA4) to the PS first crossing distribution (shown with crosses) at 2: = 0. The con'esponding absorbing baniers 
are shown in panel (fe). The crosses are overlaid on the solid curve in panel (a) because equation jA4) produces the same results as the PS case when the 
barrier lAlt is approximately equal to the scale-independent virialisation threshold, (5„ . For the barriers that dip downward, it is more likely for trajectories to 
be absorbed leftward of the corner. This depletion results in a smaller first-crossing probability rightward of the comer. 



where r — ^J2\\\/B'^ — 1 and w = yJl + 2\/B^ . The condition Q{x, y = BA) — cannot be satisfied simultaneously by both terms in 
equation l lB3t . Moreover, the se cond term does not co nverge upon applying the above condition. Discarding the second term and rewriting 
the solution in terms of r yields jMcOuinn et al .I2OO5I) 



Q{x,y) — / dr /i(r) exp 



sm{r{y-BA)}. 



(B4) 



We now apply the initial condition Q{0,y) — \B\ Soiy — yo), where 5d is the Dirac delta function, to equation iB4) to obtain ^(t) 
2\B\ e^" sin [T{yo — BA)] /tt. Integrating yields the conditional probability density 



Q{x,y\0,yo) = 



/2-KX 



exp 



B^x 



y + yo 





-iy-yof 




|exp 


2B^x 


— exp 



'{y + yo-2ABf 



2B2a 



(B5) 



The first crossing distribution can be obtained from equation l lB5b using 



Mx,AB\0,yo) = 



d 

' dS 



dS 



J^^Q{x,y\0,yo)f 
J^^Q{x,y\0,yo)^ 



ragy 



2 lay 



B <0 



where the second set of equalities follow from using equation l lB2b . For B < and B > 0, we obtain 



Mx,AB\0,yo) 



{AB_yo) 
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exp 
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